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FOREWORD 


The present report is one of a series of three reports 
describing a new computer program which predicts turbulent bound- 
ary layer behavior under conditions involving both heat transfer 
and pressure gradient. Part I serves as a summary report and 
describes the general analysis which is utilized in the numerical 
calculation scheme. In Part II the requisite low speed formula- 
tion, consisting of a constant property flow with combined 
pressure gradient and mass transfer is described. Part III 
describes the numerical and computational procedures involved 
and serves as a computer program manual. 

The titles in the series are: 

Part I - Summary Report - "Calculation of Turbulent Boundary 
Layers with Heat Transfer and Pressure Gradient 
Utilizing a Compressibility Transformation, " by 
C. Economos and J. Boccio. 

Part II - "Constant Property Turbulent Boundary Layer Flow with 
Simultaneous Mass Transfer and Pressure Gradient," 
by J. Boccio and C. Economos. 

Part III - "Computer Program Manual," by J. Schneider and J. Boccio 

This investigation was conducted for the Langley Research 
Center, National Aeronautics and Space Administration, under 
Contract No. NAS1-9624, with Mr. Kazimierz Czarnecki as the NASA 
Technical Monitor. 

The Contractor's report number is GASL TR-748. 
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SUMMARY 


This report presents the results of a theoretical investiga- 
tion of compressible turbulent boundary layer flow with heat 
transfer and pressure gradient. The analysis utilizes a com- 
pressibility transformation but differs from previous work in 
that higher order closure rules are utilized to complete the 
transformation. Specifically, by requiring that the momentum 
equations in differential form be satisfied at the wall and at 
the sublayer edge, correspondence rules are obtained which relate 
the variable property (VP) flow of interest to a constant property 
(CP) flow in which mass transfer and pressure gradient occur 
simultaneously. To implement this approach a new CP formulation 
is developed which includes both of these effects. A computer 
program based on these formulations, together with a Crocco 
integral representation for the energy field, is developed. 

Peteils of the CP analysis and the computational procedures are 
presented in companion documents. Numerical results for a variety 
of cases are presented in this report. Comparisons with earlier 
forms of the transformation and with experiment are also included. 
For the zero pressure gradient case some differences between the 
various predictions are observed. However, for the several press- 
ure gradient cases which are examined, the results obtained are 
found to be essentially identical to those given by first order 
closure rules; i.e., fc>y a form of transformation which relates 
the VP flow to a CP flow with pressure gradient but zero mass 
transfer. 
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CALCULATION OF TURBULENT BOUNDARY LAYERS WITH 

HEAT TRANSFER AND PRESSURE GRADIENT UTILIZING A 

COMPRESSIBILITY TRANSFORMATION 

PART I: SUMMARY REPORT 

By C. Economos and J. Boccio 
General Applied Science Laboratories, Inc. 


I . INTRODUCTION 


In a superficial sense, the overall objective of the current 
investigation is identical to that of previous work reported in 
Reference 1. That is, both investigations address themselves 
to the same problem area (VP turbulent boundary layers) and both 
utilize a "compressibility transformation" approach to treat this 
problem. Nevertheless, the objectives differ in a very sig- 
nificant way since the present investigation is particularly 
concerned with development of a more self-consistent form of 
transformation which eliminates some of the anomalous behavior 
observed in the earlier work. 

For a complete review of these deficiencies, the reader is 
referred to Reference 1. For the present purpose, it suffices 
to say that a major portion of the work to be described here 
involved development of a more general form of transformation 
which would allow introduction of an additional free parameter. 

As outlined in Reference 2, such a generalization would 
imply a correspondence between the general VP flow case, and a 
CP flow with mass transfer and pressure gradient such that, 
even for the zero pressure gradient VP case, both mass transfer 
and pressure gradient would persist in the CP flow. Such a 
result appears to be attractive a priori since mechanisms are 
thereby introduced which might be expected to compensate for 
the two most striking deficiencies observed in the transformation; 
namely distortion of the wake portion of the velocity profile . 
(Reference 3) and overestimation of skin friction coefficient at 
high heat transfer rates (Reference 4) . This expectation is 
based, first of all, on the demonstrated improvement in skin 
friction prediction (Reference 2) which has been achieved for 


the high heat transfer case by use of a form of transformation 
which relates constant pressure VP flow with heat transfer to 
constant pressure CP flow with mass transfer. So far as the 
wake distortion is concerned, no corresponding development has 
been reported.* Nevertheless, since the wake distortion can 
qualitatively be associated with the existence of a pressure 
gradient (c.f.. Reference 3), it may be anticipated that the 
form of transformation proposed here will simultaneously provide 
. improvement in both areas of concern. 

• * In order to exploit the proposed modification it is nec- 

'essary, of course, to have available a suitable CP formulation 
which describes turbulent boundary layer development under the 
influence of both mass transfer and pressure gradient. At the 
inception of this program such a formulation did not exist and one 
has been developed during the current study. The details of this 
development, including presentation of numerical results and com- 
parisons with experimental data, are given in Part II of this 
report. Subsequent to this development, a similar analysis was 
carried out by other investigators (Reference 5) . Comparison 
of the two methods is also included in the aforementioned com- 
panion document. 

The material presented in the subsequent sections is 
structured as follows. First, the modified closure rules which 
provide the desired correspondence are derived. We note here 
that two types are considered, one of which is rejected on the 
grounds that it does not exhibit uniformly valid behavior for 
limiting cases. Development of the working equations is then 
outlined with most of the algebraic detail either presented in 
the Appendices or omitted completely by citing appropriate 
references where such detail has previously been presented. 
Finally, a series of numerical results is presented and com- 
pared with appropriate experimental data and with predictions 
due to earlier forms of transformation. These results include 
zero and arbitrary pressure gradient cases both with and without 
heat transfer. 


* In Reference 2, an empirical correction is developed to account 
for this effect. However, this does not represent a self- 
consistent modification of the compressibility transformation 
itself. 
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Superscripts 
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III. ANALYSIS 


A. General Considerations 

Most of the detailed development of the original Coles' com- 
pressibility transformation (Reference 6) has been worked out 
many times (c.f., References 1, 7 and 8) and need not be repeated 
here. For the present purpose, it suffices to note that by 
introduction of three x-dependent stretching parameters, £,r},<y, 
for the x and y coordinates and the stream function $, a VP flow 
with fluid dynamic behavior described by 

SfiS. + Sfiv = o (1) 

dx by 


nu 




dT __ dp 
by dx 


(2) 


is transformed to a companion CP flow given by 


si + 

ax ay 


o 


pu 


au 

ax 


+ 


P v 


si = sx. _ 

ay 5y dx 


(3) 

(4) 


To obtain closure (i.e., to complete the transformation) sufficient 
restraints must be imposed to define not only the transformation 
parameters, £, r? and a, but the pressure gradient term dp/dx since 
this must be considered as an unknown in this formulation. If 
mass transfer is involved, the CP transpiration rate pv w is an 
additional unknown requiring still another restraint. Thus, for 
the most general case a total of five relations are required to 
achieve closure. 


Up until the present time this more general case has not been 
worked out in detail although in Reference 2 a possible general 
procedure was outlined. *In Coles' original work, closure for the 

* The general case has also been examined by Lewis (Ref. 9) but 
from a somewhat different point of view. A discussion of the differ- 
ences and similarities of the two approaches may be found in Ref. 2. 
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special case of zero pressure gradient and zero transpiration in 
both planes was achieved by (a) satisfying Equations (2) and 
(4) at the outer edge of the boundary layer, (b) satisfying 
integrals of these equations from the wall to the edge, and 
(c) invoking the so-called sub-structure hypothesis which 
assumes the invar ience, under the transformation, of a character- 
istic Reynolds number based on average thermodynamic properties 
within the viscous layer. Libby and Baronti (Reference 8) 
examined the pressure gradient case but with zero mass transfer. 
They also employed (a) and (b) but replaced (c) with the sub- 
layer hypothesis which evaluates the characteristic Reynolds 
number at the sublayer edge. The fourth restraint required to 
achieve closure in this case was obtained by satisfying 
Equations (2) and (4) at the wall. Economos (Reference 7) con- 
sidered the mass transfer case with zero pressure gradient. 

Closure was achieved in a manner similar to that of Reference 8 
but with the blowing rate pv w replacing dp/dx as an unknown. 

With the exception of the sub-layer (or substructure) hypo- 
thesis, all of the restraints utilized by the various investigat- 
ors can be thought of as "compatibility" conditions on the 
various unknown parameters which assure satisfaction of the 
describing equations at specified points in the flow regime. Any 
number of such conditions can be set down in a variety of ways 
so that, in principle, additional free parameters could be 
introduced with closure still being possible. For example, the 
approach used in Reference 2 to treat the case of simultaneous 
mass transfer and pressure and pressure gradient proposed utiliz- 
ing the first derivatives of Equations (2) and (4) evaluated at 
the wall to achieve closure. Alternately, Equations (2) and (4) 
could be evaluated at points other than the wall; at the sub- 
layer edge, say. 

The point of view enumerated above implies non-uniqueness of 
solution to the problem of interest. Thus, to assess the useful- 
ness of any particular choice, the entire formulation must be 
developed and implemented to the point where prediction of tur- 
bulent boundary layer behavior can be generated and compared with 
experiment. In the next section, two such possibilities are 
examined in detail. These include the use of the first derivatives 
of Equations (2) and (4) evaluated at the wall which will be 
referred to as "closure by wall compatibility." It will be shown 
that this particular choice is unsatisfactory by virtue of non— 
uniformly valid behavior for limiting cases of zero mass and heat 
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transfer and zero pressure gradient. The second choice involves 
satisfying Equations (2) and (4) at the sublayer edge. This 
will be denoted by "closure by collocation." It is found that 
this second choice yields well behaved solutions for all ranges 
of conditions. 


B. The Second Order Closure Relations 


Closure by Wall Compatibility- . - The starting point for 
deriving the closure relations is Equations (2) and (4) together 
with certain relations between the VP and CP flow parameters 
which follow from the basic transformation rules. The pertinent 
results, taken from Reference 1, are 


3_. _ p. 9_ 
9y 17 p sy 


(5) 

u =*u 
CT 


(6) 

u/u_ = uAl s u 

C 

(7) 

|5u| ^w rf j 

Su| 

(8) 

\9y/w p a ! 

t 3y 1 w 


Consider now Equation (2) evaluated at y = 0. There results 

( 5u| _ | _ dp 

^w v w \ dy/w \Sy ^ dy|w dx 

where mass transfer in the VP plane as reflected by the coefficient 
p w v w has been included for generality and it has been recognized 
that in the vicinity of the wall the shear takes on a laminar 
form. Carrying out the indicated differentiation leads to 



M ) ( — | - — 

I9y' w J \ 9y /w u 


p u 

u, r e e 

r w 


du 
€ 
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(9) 


where we have used the Euler relation oe u edu e /dx = - dp/dx. If 
the same procedure is applied to Equation (4) there results 
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( 10 ) 



which is identical in form to Equation (9) except that the 
viscosity derivative vanishes since, for CP flow, p is a constant. 
We now note that the relations (5) , (6) , and (8) imply 


B 8 U _ 0_ 

, 9y s jw if 




Bo 

5y|w 


3u 

Byjw 
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( 11 ) 


which permits elimination of the second derivatives appearing 
in Equations (9) and (10) . There results 
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where we have recognized that the wall shear is related to the 
velocity field by 


- ,3u. 

T w ■ “ ( i?>w 


, 3u , 

T w ^w By w 


In the further development it will be assumed that the thermo- 
dynamics can be related to the velocity field by means of a Crocco 
integral. Accordingly, we can write 


(&U) = (hU. 

Vw 


(*■**■) (— ) 
3u w By w 
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By w 3u w By w 


Then Equation (12) can be written 
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where the following definitions have been introduced 


Z = 


P = 


2 
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o o v 
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c_ pu 
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C f °e U e 
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and the remaining parameters U , c f , p, etc., have been defined 
in the List of Symbols. Note fhatthe terms appearing in 
Equation (14) are of three types; terms proportional to an in- 
jection rate (Z, Z) , terms proportional to velocity gradients 
(P, P) , and a term Gr which depends solely upon the thermodynamics 
of the VP fluid. W 

It is useful at this point to show how special forms of 
Equation (14) have been utilized by other investigators to obtain 
closure. In addition, the effect of the fluid thermodynamics on 
the CP flow behavior will also be examined qualitatively. The 
latter is best accomplished by considering a specific thermo- 
dynamic system in order to evaluate explicitly the parameter G i w * 

For this purpose we assume that the CP fluid is a perfect gas 
with an exponential viscosity- temperature relation and Unity Prandtl 
number. It is further assumed, for the case of injection, that the 
injectant gas is identical to the external fluid. Then the Crocco 
integral for total enthalpy implies 


w 


= (1-a) (1-W) 
W 


(15) 


where a denotes the viscosity exponent (for most gases a < 1) and 

W is the temperature ratio T w /T t . Evidently, for the adiabatic 

wall case (W = 1), Equation (15) e implies that Gi w = 0 while for heat 

transfer to the surface (W<1) , >0. Note however, for the 

special case a - 1 (i.e., a linear viscosity temperature relation) 

Gi =0. 

- L w 
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Now consider the special case Z = Z = 0, which is the flow 
configuration examined in References 1 and 8. The closure 
relation becomes, from Equation (14) 


P = P 



( 16 ) 


which implies some interesting possibilities insofar as the CP 
pressure gradient is concerned when a ^ 1. It is apparent, for 
example, that P does not vanish when P = 0 except for the adia- 
batic wall case. It is also possible, of course, to achieve the 
case P s 0 by considering a VP flow where P - =0? whether 

this corresponds to a physically realizable configuration is not 
known at this time and will not be considered further. 

The most relevant possibility, for the present purpose, is 
the situation corresponding to VP flow with heat transfer and 
zero pressure gradient. This is of particular interest here in 
view of the fact that much of the anomalous behavior exhibited 
by earlier forms of the transformation has been associated with 
experimental data obtained under such conditions. 

From Equation (16) with P = 0 we note that 

P = < 0 when G, >0. 

x w 

That is, with heat transfer to the wetted surface (G;l > 0) , this 
particular_closure rule implies the existence of adverse pressure 
gradient (P < 0) in the CP plane, while for the "hot wall" case 
(G]^ < 0) a favorable pressure gradient occurs. The first of these 
results appears to be qualitatively inconsistent with observation, 
since, as pointed out in Reference 3, the distortion of the wake 
component is such as to imply a favorable pressure gradient. Note 
however that this distortion has been observed even for adiabatic 
flow for which this particular formulation yields no effect. 

A somewhat parallel situation occurs for flow with mass 
transfer but zero pressure gradient as treated in Reference 7. 

The closure relation in this case was taken to be 

Z = Z + Gx (17) 

J -w 

which follows from Equation (14) when P = P = 0. It is now easy 
to see how the use of two different closure rules yields different 
solutions for the same physical problem. That is, the zero mass 
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transfer, zero pressure gradient, VP case can be associated either 
to a CP flow with pressure gradient given by 

P = - Gi (18) 

x w 

or with constant pressure flow with transpiration (or suction) 
given by 


Z 



(19) 


The second of these formulations, of course, contributes nothing 
in the way of explaining the wake distortion. However, as 
demonstrated in Reference 2, it can improve skin friction pre- 
diction for the cold wall case G;l >0. It might be anticipated 
that some intermediate combination of these effects would provide 
improvement in both areas, which is, of course, the motivation 
for the current investigation. 

Returning now to the general problem, we require an additional 
relation involving the parameters P and Z. Here we utilize the 
"second wall compatibility" relation obtained by differentiating 
Equations (2) and (4) with respect to y and y, respectively, 
and evaluating these at y, y = 0. 

This procedure yields 
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the third derivatives: 
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Now combine these last three equations to obtain 
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where we have also utilized (7) , (8) , and (13) . To proceed, it 

is necessary to relate the last term on the right-hand-side of 
Equation (20) to the velocity field. Using Equation (13) it is 
easy to show that 
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The second derivatives can then be eliminated as before yielding 
the desired result which is 


z(z-p) = (z+g 2 ) (z+g 3 ) + g 4 - (z+g 3 ) P 


where 
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( 21 ) 


Note that all of the G^ 1 s depend only on fluid thermodynamics. As 
with the term G lw , they can also be expressed in terms of 
specified flow parameters according to 


14 



1-W 


( 22 ) 


G„ - - 


G 


3 


G 


4 


(4-3a) 


1-W 

W 


( 1 +a-a 2 ) 


(irW ) 2 
v w ' 


2 ( 1 - rc .) ( V -» Ma / 2 

W 1+(v-1)M 2 /2 

e 


(23) 

(24) 


It is interesting to note the Mach number dependence of G 4 and 
the fact that it does not necessarily vanish for the adiabatic 
case . 


Equations (14) and (21) represent the two additional re- 
straints required to treat the general VP problem involving both 
mass transfer and pressure gradient. They should, however, be 
equally applicable and uniformly valid for limiting cases 
(i.e., Z -► 0, P -♦ 0, W -» 1) . In order to examine this limiting 
behavior it is convenient to introduce some new notation. Let 
Z = 0; then Equation (21) reduces to 

Z(5-P) = G-GP (25) 

5 3 

where 

°5 2 G 2 G 3 + G 4 (26) 

Several limiting cases can now be considered as itemized in 
Table I. The first entry therein corresponds to the zero pressure 
gradient adiabatic wall case. From Equations (15) and (23) it 
follows that G^ = G 3 = 0 as has been indicated. Equations (24) 
and (26) imply ¥hat G 5 ^ 0 (assuming M e ^ 0) . In this circumstance 
Equations (14) and (25) imply 

Z - P = 0 

Z(Z-P) ^ 0 

The only way in which these equations can be satisfied simultan- 
eously is for Z and P both to be infinitely large which is clearly 
unacceptable from both a physical and computational point of view. 
This irregular behavior could be suppressed by taking the viscosity 
exponent to be unity. The resulting situation is indicated by the 
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second entry in Table I. In this case we find that Z = P but that 
neither one is explicitly defined. Even if we consider this an 
acceptable situation, the irregular behavior has not really been 
eliminated. This is demonstrated by the third entry wherein 
the adiabatic case with non-zero pressure gradient is considered. 

Again, it is found that in_order to satisfy the closure relations 
(14) and (25) , both Z and P must take on infinitely large values. 

In view of these results it is concluded that a consistent set of 
.closure rules cannot be obtained in this manner, at least for the form 
of transformation utilized here (e.g., Eqs . (5) , (6) and (8) . Accordingly, 

an alternate approach is required and this is discussed in the next 
section. 

Closure by Collocation . - In this approach Equation (14) 
is also utilized. However, Equation (21) is replaced by satisfy- 
ing Equations (2) and (4) at the sublayer edge. This would yield 
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where subscript s indicates that the variables have been evaluated 
at their respective sublayer edge. In order to be consistent with 
the CP formulation, however, the following approximations are 
appropriate : 


d£ 

dx 


P U C 


<f§) 

3x s 


« 


(ki) 

3y ; 


- pv 


v ay s 


V ^ V 

s w 

By analogy we assume that similar approximations are appropriate 
for the VP flow. Thus if it is further assumed that p s P w and 
that at the sublayer edge the shear is due solelv to laminar viscosity. 
Equations (27) and (28) reduce to 


( 


B 3 u . 
By 2 s 


rVw 

1 ^s 


— (**) 1 (^) 

p s V's iky's 


, B a u x _ pV w 

'ay 8 * s jj 


By s 
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By utilizing relations similar to Equation (11) and (13) and 
eliminating the second derivatives between these last two 
equations, there is obtained 


Z 




<31 


s 


where 


G ls 


, 5-tnpg 
y bu ’s 


and 


- _ -,du. 

T s * 


,au. 

" ^s dy s 


For the further development it is convenient to express the 
laminar shear r s in terms of other variables of this formulation. 
For this purpose we can use Equation (10) of Part II of this 
report to show that 


/ 

T / T 
S W 



r7 

t-X 


where u denotes the velocity ratio occurring at the edge of the 
sublayer. We will also restrict further consideration to the 
case of zero mass transfer in the VP plane. Accordingly, the 
desired closure relations become 


Z - P = G. - P 
lw 

1+u 2 G 3 -s 
s 

where, to be consistent with Equation (15) , G-^ would be 
expressed as 5 

(l-W)+2m e U g 

G 1 = ( 1-a ) w+(l-W)T3 +1ft H 2 
s s e s 


(29) 

(30) 


(31) 
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where 


(y — 1 ) M 2 / 2 

~ _ G 

m e l+(y-l)MV2 


Qualitatively, Equations (29) and (30) appear to be well behaved 
for all conditions of interest. We note especially_that even 
for zero pressure gradient adiabatic flows neither Z nor P 
vanish which would certainly be the type of effect being sought 
here. Note also that as M - 0 and W - 1, there is obtained 
Z = 0, P = P which impliel that the system reduces uniformly to 
the CP problem as it should. 


A more quantitative assessment of the general behavior of 
these relations is not possible without full implementation of 
the entire formulation. This is due_to the fact that the velocity 
ratio u is a function, not only of Z, but of the CP skin fric- 
tion colfficient as well. Nevertheless, some representative 
results have been generated for the zero pressure gradient case 
by considering u as a parameter. That is, we let W, M , u g and 
a vary over reasonable range of interest and evaluate tBe 
corresponding values of Z from 

G l c 


Z = 


1-u G-i 
s • L s 


(32) 


which follows directly from Equation (30) . Some representative 
results of this procedure are shown in Figures 1 and 2 . Figure 1 
shows the effect of Mach number and wall temperature ratio at a 
fixed value of a with the velocity ratio u as a parameter. In 
Figure 2, u is maintained constant with a taken as the parameter. 
As may be noted, Equation (32) is well behaved throughout a wide 
range of all of the parameters even up to and including infinitely 
large Mach numbers and vanishingly small wall temperatures. By 
contrast, the first order closure rule represented by Equation 
(19) , which incidentally does not include any Mach number depend- 
ence, exhibits singular behavior when W -> 0 as may be seen by 
examining Equation (15) . 

Superficially, it would appear that the second order closure 
rules derived here represent a means for resolving all of the 
difficulties previously encountered with the transformation. 

Closer examination of the data shown in Figures 1 and 2 implies 
however that this may not be the case. This is particularly 
evident by noting that for the adiabatic wall case W = 1, negative 
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values of Z occur for all non-zero Mach numbers. But from 
Equation (29) with Gi = P = 0, we have Z = P < 0 which implies 
the occurrence of adverse pressure gradient in the transformed 
plane. Since this is precisely opposite to the effect desired, 
the likelihood of improving the performance of the transforma- 
tion approach by second order closure is substantially dimin- 
ished. This will be borne out by the results presented in 
Section III- 


C. Working Equations 

General Remarks . - The final result of this section will be 
a system of ordinary differential equations whose matrix repre- 
sentation can be written 


Al 

A9 

A17 

A25 

A3 3 

0 

0 

- 

0 


an R) ' 


0 

A2 

A10 

A18 

A26 

A3 4 

0 

0 

0 


an (p) ' 


C2 

A3 

All 

A19 

All 

A3 5 

0 

0 

0 


(ir) 7 


C3 

0 

0 

0 

A28 

U 

0 

0 

u 


{ » ~ TT^ ' 

= 

C4 

0 

A13 

0 

A29 

A3 7 

0 

0 

A63 


(V*/)' 


0 

0 

A14 

0 

A30 

A3 8 

A46 

0 

A65 


an a) ' 


0 

0 

0 

0 

A31 

0 

A47 

A55 

A62 

i 

| 

an n) ' 


0 

0 

0 

0 

0 

0 

A48 

A56 

A64 

j 

ix>' j 


C8 


It is believed that by setting down Equation (33) initially, the 
detailed development which follows will be facilitated and 
clarified. This procedure will also provide an opportunity to 
compare the present formulation with that of Reference 1, at the 
outset. By virtue of this comparison, those similarities which 
occur in the two methods can be exploited to reduce much of the 
tedious algebra involved in the derivation of Equation (33) . 

In addition to this derivation, appropriate initial 
conditions will be developed in this section. For this purpose 
it will be necessary to outline a procedure for deducing skin 
friction coefficients from experimental velocity profile data 
according to the second order closure rules employed here. 
Similar procedures for the various other closure rules^will also 
be set down. Although these so called "Clauser Plots" have been 
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utilized by previous investigators (c.f. References 3, 7, 10) 
they have not, to the authors' knowledge, been set down in a sys- 
tematic manner in any of the literature. 

Finally, we note that the reason conventional notation has 
not been used in the coefficient matrix of Equation (33) is 
twofold. First of all it was desirable to avoid any possibility 
of confusing the various terms appearing here and in Reference 1. 
Secondly, the choice of notation anticipates the logical flow 
of certain "canned" subroutines used by the CDC 6600 computer 
for which this analysis has been programmed. A further discussion 
of this aspect may be found in Part III of this report. 

Comparison of First and Second Order Closure Formulations 
for the Pressure Gradient Case . - In Section HI .B, it has been 
pointed out that closure was achieved in Reference 1 by means of 
Equation (16) . The resulting formulation has the matrix 


representation 

(c.f. 

Equation 

(56) of 

Reference 1) : 


— 








""" — 



An 

a 12 

A 13 

A 14 

0 

0 

0 


(tn 0 ) ' 


0 

A 21 

A 22 

A 23 

A 24 

0 

0 

0 


( iir)' 


C 2 

A 31 

A 32 

A 33 

A 34 

0 

0 

0 


Un R)' 


C 3 

0 

0 

0 

A 44 

0 

0 

0 


(In Ue) ' 


°4 

A 51 

0 

0 

0 

A 55 

0 

A 57 


(tn a) ' 


0 

0 

0 

0 

^4 

^5 

*66 

A 67 | 


(in 7/) ' 


0 

_0 

0 

0 

0 

A 75 

A 76 

^11 - 


.(x)' 


_ C 7 _ 


The most obvious difference between Equations (33) and (34) 
is that the former is an eight-dimensional system as compared to 
a 7 x 7 for the latter. Evidently this increase in dimension 
is due to the introduction of the additional free parameter 
corresponding to transpiration in the CP plane. In Eq. (33) 
this parameter is represented by the dependent variable v w whose 
relation to the parameter Z previously introduced will be 
indicated in the subsequent development. Aside from this rather 
essential difference the two systems bear a considerable similar- 
ity in the sense that the same fundamental relations are utilized 
to derive the final working system of equations. For Equation 
(34) the specific relations, proceeding in order from the first 
row, were: 
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(a) the CP skin friction law (38) 

(b) the CP momentum integral equation (34) 

(c) the CP auxiliary equation (43) 

(d) the first order closure relation (26) 

(e) the sub-layer hypothesis (45) 

(f) the edge compatibility condition (21) 

(g) the VP momentum integral equation (30) 

where the numbers in parenthesis refer to the appropriate equa- 
tions as given in Reference 1. Note that item (d) is identical 
to Eq. (16) of this report. 

The first point to be made regarding Equation (34) is that 
the first three relations (a) through (c) constitute the describ- 
ing equations for the CP flow. If the CP pressure gradient 
parameter* (£nUe) * were a known function, these three equations 
could be solved independently of the remaining four equations. 

Here, the two groups of relations are coupled by virtue of the 
fact that (£nU») must be considered as an unknown. As for this 
second group oi relations, they represent the four restraints 
which are required to define the three stretching parameters of 
the transformation and the unknown CP pressure distribution. 

The analogy with the present formulation should now be 
apparent. Again, a group of relations is required to describe 
the CP flow, but, of course, in this case they would involve 
both a pressure gradient and a mass transfer parameter. In 
Part II of this report such a system is developed and it is found 
that, once again, three equations are involved which are com- 
pletely equivalent to items (a) , (b) and (c) above. The second 

group of equations must now involve five rather than four, funda- 
mental relations to account for the additional unknown parameter 
V . Of these, three are exactly equivalent to items (e) , (f) and 

(g) . The remaining two are provided by 1 the closure rules developed 
in the previous sections; i.e., Eqs. (29) and (30). 


* U e bears the same relation to the parameter P that V w does to 
Z. The explicit correspondence between these variables will 
be indicated later. 
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With this as a background we are now in a position to 
proceed with the derivation of Equation (33) in a straight- 
forward manner. This will be accomplished in the next three 
subsections. 

In the first of these the working form of the two closure 
relations, which are unique to this formulation, will be derived. 
Then the relations corresponding to items (e) , (f) and (g) will 

be considered. Finally, the CP formulation of Part II of this 
report will be arranged in a manner suitable for numerical 
integration and incorporation into the matrix representation of 
Equation (33) . 

Working Form of the Closure Relations . - The working depend- 
ent variables which are employed in the numerical integration 
scheme have been indicated in the solution vector appearing in 
Equation (33) . The first five members of this vector all 
represent CP variables and have been defined in Part II of this 
report. Furthermore, with the exception of the mass transfer 
parameter V w , all of the variables are identical in meaning to 
those utilized in Reference 1. Accordingly, to proceed with this 
derivation, it is only necessary to show the relation between Z 
and P and these working variables. From the definition of these 
latter and those appearing on page 11 of this report, it follows 

_* 

P = UnU e ) ' 
e 

can be written 


or, taking into account the definitions of P and the correspond- 
ence* between the CP and VP skin friction coefficients 

V U G-i 

, . - , , w e x w , 

UnUe) = , -r- r~ + 

<p (P 


V r 




d-lnU 
£ 

d x 


(35) 


immediately that 
5 2 V 

5 -V s ’ 

e 

Accordingly, Equation (29) 

V 

■=§■ - UnU e ) 
_ <P 



u 

e 


* Correspondence between the CP and VP variables as implied by 
the transformation are listed in Appendix A. 
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This equation corresponds to the fourth row in Eq. (33) with 
A28 equal to unity and C4 taken equal to the right hand side 
of Eq. (35).* Complete evaluation of C4, of course, requires 
specification of both the external boundary conditions ^ ^ 

(d-tn U /d\) as well as the wall boundary conditions (G ^ W »P W »4 W ) 
the latter also involving the thermodynamic behavior 
of the VP fluid. Detailed evaluation of these parameters is 
discussed later in this section. 

Consider now the second closure rule, Eq. (30). Proceeding 
as before, we can write 

tp 2 v 

H e [ 1 - (1 + u ) G, ] = 0 (36) 

1 U e s Is 

Since this is an algebraic equation it is necessary to differ- 
entiate it for the purpose of incorporating it in Eq.(33). In 
carrying out this differentiation the implicit dependence of 
u_ and G., on the dependent variables must be recognized. From 
Eq. (20at of Part II of this report it can be shown that 

u = — —■ [ exp (v+y*) - 1] (3 7) 

s — + w s 

<pv 
^ w 


where v = « — 
w U 

e 

and y + = fcn(v + ) as determined from the transcendental 
equation (20c) of Part II. Accordingly, 

u = fcn($,U ,V ) 
s e vr 

In view of Eq. (31) we may also anticipate that in the 

general case, G will exhibit a dependence on both u and the 
■“■S s 

* All of the coefficients appearing in Eq. (33) are defined 
explicitly in Appendix B. 
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streamwise coordinate X through a corresponding Mach number 
variation. Eq. (36) is then of the form 

H (<p U , V , X) = 0 
1 ' e w 

Thus, differentiation of Eq. (36) yields 

bR a-H dH SH 

(-_!_) {ln(p) . + ( 1-) Un u ). + (— i) (y ). + (T~ i )(x)' = o 

a-tnS a-tnu e av w x 

e w ( 36a ) 


which corresponds to the fifth row of Eq. (33) where the 
coefficients A13, A29, A37 and A61 are related to the partial 
derivatives in an obvious way. The latter are evaluated for 
a particular choice of thermodynamic behavior in Appendix C. 

Working Form of the Remaining Transformation Restraints : - 

(a) Sublayer Hypothesis - The basic form of this restraint is 
given by 


p u y 
s s s 


pu y 
s s 


(38) 


By utilizing the various correspondences between the CP and VP 
variables given in Appendix A and the definitions this can be 
written - + 



s 1 


~ — 
"s y s 



(39) 


The dependence of this relation on the thermodynamic behavior 
of the VP fluid is apparent. Proceeding as in the case of 
Eq. (36) we recognize that Eq. (39) is of the form 

H (^, U ,V ,o,x) = 0 

2 e w 
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Differentiations yields 


94nH 2 

b in,/) 


__ ^ 5-C,nH 2 

(t ”* ) ‘ + 8^0" 
e 


(tnU ) ' + 
e 


d^nH 2 


bV 


w 


~ bl n H 2 

(V ) ' t (ina) ' + — x'=0 

w 


(39a) 


which corresponds to the sixth row of Eq. (33) . Explicit eval- 
uation of these partial derivatives is presented in Appendix C. 


(b) 

Eq. 


Edge Compatibility - This restraint follows directly from 
(6) evaluated at y, y -* ®. There results 


e a e 


or, in terms of the working variables 


U 

e 


^rr 

a 




u 

e 


f a n\ 

f 


The working form of this relation, corresponding to the seventh 
row in Eq. (33) follows by differentiation. There results 


dfcnU /v 

Unuj •- Una) ' + Unrj) * — ^ — ~ x’ = 0 


(40a) 


where we have recognized that w and (a/ri) Q are constants. 

o 

(c) VP Momentum Integral - Since this formulation will corre- 
spond to the case of zero mass transfer in the VP plane, the 
momentum integral equation remains identical to that used in 
Reference 1 and can be written 


d£ 

dx 


+ 



u dx 
e 


(H+2) 


+ 



o dx 
^e 



(41) 


The CP counterpart of this equation must, of course, reflect the 
existence of mass transfer. Accordingly, it takes the form 
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(42) 


dg 

dx 
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u dx 
e 


(H+2) 



+ 


v 


w 


U 

e 


After introduction of the various definitions these can be 
written 


dtn0 + 

dX 

(tn0) ' 


(H+2) 


dtnU 
e 

dX 


dtno 
€ 

dx 


+ (H+2) (tnu ) ) = 5"('=s' 
e 0 <p 


c u -1 

— (— ) 

26 v o 

e 


V 


u- 


+ (— ) 

u e * 


-1 


(41a) 


(42a) 


Now the correspondence between 0 and 0 according to Appendix A 
is given by 0 = 0 /t). Accordingly 

(tn0 ) = -, [(tn0) '-Unri) 1 1 

d X X X 


Hence, Eq. (41a) can be written after some rearrangement 


(H+2) (tnuj ’ + (tn^) '- x '{(H+2) 
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* + 


dX _ d x 
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e . _f ( _e } 1 = 

26 V o J 


V e -1 
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e V V * 


It is also convenient to eliminate (tnU ) 1 by means of Eq. (40a) . 
There results e 


dtnU 


dtnp 


(H+2) Una) ’-(H+l) (tnr?) '-{(H-H)- ^ 6 -(H+l) 6 +(H+2) 


d^njLt 


d X 


c u -1 
20 


u -1 
e 


u -1 -j -j V e 

/ e. \ , 1 /l . w w o. 

(“) X r(=F + jT~) (—=-) 

V O J A 0 (0 U v 


(43) 


In Part II of this report, a CP Reynolds number based_on boundary 
layer thickness, 6, is defined according to R = u e 5/v which is 
one of_the working variables appearing in Eq. (33?. Let the 
ratio 0/6 be denoted by £. Then 
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1^ 
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-1 
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20 V o 
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PM cr 
'w w 
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where we have utilized several of the correspondences given in 
Appendix A. Note that by introducing the parameter R the 
dependence of this system on the unit Reynolds number of the CP 
flow and, consequently, on its thermodynamic state, has totally 
been eliminated. 


Equation (43) can now be written 


d£nU _ d£np 

(H+2) UnZ) 1 - (H+l) (■tnri) { (H-H) — e -(H+l) — + 

L d x d x 


(H+2) 



<p s ER 




(44) 


In Eq. (44) the \ derivatives of U , p , and /i , are, of course, 
considered to be specified and hence are known® Furthermore, 
the CP formulation provides relations of the form 


L = Z (<p, ir, U , V ) 
e w 


H = H (<p, n, u , V ) 
e w 


Consequently, provided only that the VP form factor H is ex- 
pressed in terms of the working variables, Eq. (44) is in a form 
suitable for incorporation in Eq. (33) . This representation 
must await specification of a particular thermodynamic model for 
the VP fluid since, in accordance with the correspondence (A8) 


H 


H + 




(i -1) dy 
P 


(45) 


The integral appearing in this expression is evaluated later in 
this section. Equation (44) is represented by the last row of 
Equation (33) . 
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Working Form of the CP Formulation . - The describing 
equations for the CP formulation are given by Eq. (38) of Part II 
of this report. It takes the form 


A 1 
A 2 
A 3 


A 4 
* 5 
A 6 


A 7 
A 8 
A 9 


(•livo) ' 
UnR) ' 
(ir) 1 



(46) 


where derivatives of V and U are included in the forcing 
terms C.. For the CP problem 6 itself this is appropriate since 
these are considered to be given forcing functions. For the 
present application they represent unknowns so that this system 
has to be modified accordingly. The procedure is straightforward 
and involves only the transposition of the coefficients of (V ) ' 
and (^.nU ) ' from the right hand side of Eq. (46) to the left 
hand side. After some rearrangement there is obtained the first 
three rows of Eq. (33) where the coefficients in (33) and (46) 
are related by 


Al = A 4 , A2 = A g, A3 — A 5 

A9 = A lt A10 = A 3 , All = A 2 

A17 = A 7 , A18 = A g, A19 = A 8 

A25 (l,nU ) ' + A33 (V ) ' = - C 
e w x 

A26 (-InU ) ' + A3 4 (V )' - C2 = - C 3 
e w 

A27 (in U ) 1 + A35 (V )' - C3 = - C 2 
e w 


For convenience in the further development it is useful to 
note here that the first of these equations arises by differentia- 
tion of the algebraic equation corresponding to the CP skin 
friction law; (Eq. (28b) of Part II) . This equation can be written 


(1 + <p s 
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Thermodynamic Representation . - The influence of the VP 
fluid thermodynamics on Eq. (33) is reflected by several para- 
meters appearing in Eqs. (33), (36), (39), and (45). These 

involve viscosity and density ratios at the wall and sublayer 
edge (p^, p g , ^ g ) and derivatives and integrals thereof, i.e.. 
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1 
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' 9(1 ' w 


= (temji) 

1 1 9u 

s 



_ + 
dy_ 

P 



Each of these must be related to the velocity field and sub- 
sequently to the working variables in uiuei' to complete the 
formulation. As indicated in Section HI .B, a Crocco integral 
representation will be utilized for this purpose. For a per- 
fect gas, taking into account that the pressure across the boundary 
layer is uniform, this implies that 


T 


1 

P 


W + (1-W) u + m u 2 
e 

1 + m 

e 


(48) 


It will also be assumed here, for additional generality, that 
the VP fluid viscosity can be represented by a Sutherland type 
law 
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3/2 


T+198.6 


Then it can be shown that 
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(49) 
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(50) 
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(52) 


(53) 


Note that since T = T^_ (1+m ), Eqs . (51), (52) and (53) require 

specification of the externa? stream stagnation temperature, which 
is, of course, taken as constant in this formulation. 


To evaluate the first density integral we make use of the 
velocity distribution given by Eq. (20a) Part II. Taking into 
account Eq. (48) there results 




W+(l-W) (X 1 / ( p) 

— fV 

1+m 
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m (X 0 V 

e 2 ^ 


(54) 


where 



are functions of v w + only and are defined explicitly in 
Appendix D. 


Evaluation of the second integral is somewhat more straight- 
forward since the final result can be expressed in terms of pre- 
viously defined parameters. That is, using Eq. (48) it can be 
shown that 
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(55) 
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Thus, from Eq, (45) the correspondence for the form factor 

becomes _ 

WH-m 

H = - 
1+m 

e 


( 56 ) 


Boundary Conditions . - (a) Edge Conditions - Inviscid 

Flow: So far as the computer program is concerned, the edge 

conditions are imposed by inputting a table of Mach number versus 
streamwise length normalized with respect to a reference length 
l . The internal logic of the program then smooths these data, 
if desired, and also numerically evaluates the first derivatives 
thereof. The details of this processing is given in Part III 
of this report. Here, we will indicate only how the various 
inviscid parameters appearing in Eq. (33) are related to this 
Mach number table. 

For a calorically and thermally perfect gas the isentropic 
relations imply that 


~ h 

u = M (1+m ) 
e e e 


(y RT t ) 


(57) 


Accordingly, since the stagnation speed of sound is considered 
to be constant here 


or 
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(58) 


where M and x/Z are the tabulated input values and Re denotes 
the initial unit Reynolds number (u / v ) • Note that m deriv- 
ing Eq. (58) we have made use of th§ rllation 


dm 
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d X 


2m (1+m ) dM 
e e e 


M 


dx/z v Re l ‘ 
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which follows directly from the definition of m 


(59) 
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Proceeding in a similar fashion it can be shown that 
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In view of Eq. (60) a Sutherland-type viscosity law would then 
imply 


d X 


2m dM 
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M dx/t 1 2 
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e to 

e 


(62) 


The derivative of the kinematic viscosity v follows directly 
from Eq. (61) and (62) since 6 


d.t n ty e d£n^ e d-tnp^ 

d X d X " d X 


(63) 


Note that several of the foregoing relations are also required 
to evaluate the partial derivatives dH^/9x and S-t-nH^/dx appearing 
in Eqs . (36a)and (39a) respectively. 

(b) Wall Conditions: A detailed examination of this system 

reveals that the state of flow at the wall is completely 
characterized by the parameter W which represents the ratio of 
wall temperature to external total temperature. This parameter 
is imposed in a manner identical to that used to input the ex- 
ternal Mach number distribution. That is, a table of W vs x/£ is 
read in, smoothed and differentiated as discussed in Part III. 
Thus for the present purpose it is only necessary to recognize 
that 


dW _ dW # _1 

dx <3cx/l Re o^ 


(64) 


Initial Conditions . - A detailed examination of the system 
represented by Eq. (33) reveals that 


32 



a) The initial value of the independent variable 
X is arbitrary. 

b) The initial value of x is arbitrary except that 
the choice must be consistent with the input 
Mach number and wall temperature tables. 

c) The initial value of U and 7f are arbitrary, 
as is rj Q , but they must be selected in a 
consistent manner; i.e., they must satisfy the 
algebraic equation (40) ^ This is most conveniently 
accomplished by taking U = rj =1. 

o 

The five remaining dependent variables R, ( o, tt, V and a, are 
connected by three algebraic equations corresponding to the CP 
skin friction law, Eq. (47), and Eqs. (36) and (3 9). Accordingly, 
only two initial values need to be prescribed.* These will be 
taken to be <p and jr . The reason for this choice will be made 
clear by virtue of the ensuing development. 

As in Reference 1, two types of initial conditions will be 
considered; namely, conditions corresponding to a "leading edge" 
and conditions associated with some downstream station where 
boundary layer profiles are prescribed. 

Since the formulation developed here corresponds to fully 
developed turbulent flow, initialization at a leading edge is 
necessarily an approximate concept. Nevertheless, if the concept 
is to be employed at all, it should at least exhibit qualitatively 
consistent behavior in the sense that the solutions sufficiently 
far downstream should be independent of the particular initial 
values selected. 


* As far as the computer program itself is concerned, it will 
accept arbitrary initial values for all of the variables in 
order to permit restart and continuation of any particular 
calculation. Such a procedure is admissible provided a self- 
consistent set of values is utilized. If this is not the case, 
the inconsistency will show up by means of certain output 
parameters which provide a running check on whether the alge- 
braic equations are being satisfied as discussed in Part III 
of this report. 
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To establish appropriate initial values in line with this 
point of view, a series of numerical calculations were made. Some 
representative examples are shown in Figure 3. Here the varia- 
tion of skin friction with streamwise distance corresponding to 
various choices of initial conditions are shown. It is apparent 
that the solutions exhibit virtually no dependence on the initial 
choice of v for values in the range associated with zero pressure 
gradient CP flows. On the other hand, although all solutions 
eventually merge, some effect of the initial choice of <r> persists 
in the near region. Based on the criterion that the solutions be 
identical for \ & 10 , a suitable choice of <p appears to be 
( p ) =15. Thus, the "leading edge" option which has been provided 

in the computer program utilizes the values 7 t q = 0.6, cp Q = 15. 

We now wish to consider initialization at any downstream 
station where profile data have been prescribed. This requires 
the use of so called "Clauser Plots" which were mentioned 
briefly in a previous section. This procedure was first used by 
Clauser (Reference 11) in conjunction with analysis of CP flow 
profiles obtained on impermeable surfaces. For this case, the 
procedure is a relatively simple one and it is worthwhile to 
review this elementary case briefly here. 


We suppose that the velocity distribution at some stream- 
wise station has been determined experimentally in the form 
u vs R_ where R_ e u y/i>. From this profile it is desired to 
infer £he assoc^ated e value of skin friction coefficient c ^/2 = l/(p 2 . 
This is accomplished by comparing a semi-log plot of the experi- 
mental profile with a series of curves generated from the equation 


u 


1 

<p 




(R-/c 0 ) 


(65) 


using several arbitrary values of the parameter Equation (65) , 

of course, is the so-called "law of the wall” corresponding to 
CP flow without mas| transfer and is obtained from Eq. (20b) of 
Part II by taking v = 0. From this comparison the value of <n 
which best fits the given data can be determined by inspection.. 

An example of such a "Clauser Plot” is shown in Figure 4. Note 
that the shear is quite unambiguously determined by this method 
and that it is in good agreement with the value obtained directly 
by means of a skin friction balance. 


34 



Generalization to the CP case with mass transfer is due to 
Stevenson (Reference 12) and has been used by others (References 5, 
7) . In this case the velocity distribution is compared with 
curves corresponding to Eq. (20b) , Part II, which can be written: 


u = - { 

<p l 


k 2 (1+k 2 V> 


i^ !1+ 2 1 tn( y*’ 


? 


+ 4^; V < V*’ } (66) 

Since V is known in the case of CP flow, Eq. (66)^is exactly 
equivalent to Eq. (65) in the sense that a_unique u vs R- varia- 
tion results for any particular choice of 0 . A representative 
example of such a Clauser Plot has been given in Figure 1 of 
Part II and is repeated here as Figure 5, for convenience. Again, 
the determination of <p (and therefore of c^) is clear cut. 

Once cn has been determined, it is also possible to assess 
the value of the profile parameter ir for any particular profile. 
For this purpose the relations 


IT 





(67) 



(1+V 0 a ) 
w 


% 



( 68 ) 


are utilized for the zero and non-zero mass transfer cases, 
respectively. These relations follow from the skin friction law, 
Eq. (47) , where we have recognized that 


RU s R- b u b/v 
e 6 e 


which can also be determined from the experimental profiles. 


* Note that U e has b een take n equal- to 

appropriate since the ultimate purpose of this procedure 
is to establish initial conditions. 
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The equivalent method for impermeable VP flow was developed 
by Baronti and Libby (Reference 3) . In this case, however, the 
procedure is somewhat more involved since the measured profiles 
in the form u vs R = u y/u must first be transformed to the 
corresponding CP f&rm. e Since u can be interpreted as either the 
VP or CP velocity ratio (c.f. Eq. (7)) no stretching of the 
ordinate is required. On the other hand, points corresponding 
to various values of R must be transformed according to 


R- 

y 







d R 

y 


d R 

y 


(69) 


as indicated in Appendix A. Thus, if the density variation has 
also been obtained experimentally, a semi-log plot of u vs R-/cr 
can be prepared from the measured data. This must then be com- 
pared with appropriate forms of the "law of the wall." 


Depending upon the type of closure utilized, the u vs R- 
variation would be obtained from either Eq. (65) or ( 66 ) . ^ 

Note however that the abscissa must be divided by a corresponding 
value of a before the comparison can be made. Furthermore, if 
Eq. ( 66 ) is utilized, a consistent value of V must also be 
selected. The details of how this is accomplished for each of 
the various closure rules will now be outlined. 


"Zero Order" Closure (ZC) : This form of transformation is 
utilized in Reference 3 and presupposes that V = dp/dx = 0. 

In this case, the appropriate value of cr which is required would 
be obtained from Eqs. (39) and (54) specialized to the case V =0. 
There results ' 


p W+ (1-W) (5.3/<p) + m (3 7 .45/co s ) 

S 0 


u 1 + m 

*s e 


(70) 


where p g and /j are evaluated from Eq. (50) and (51) with 
u = 10. 6 /( 0 .* For prescribed values of m and W, Eq._(70) is 
ol the form a = cr (<p) . Accordingly, for any choice of ip, the 

* Tn Appendix E, it_is shown that as V - 0, Eq. (3 7) implies 

u = y s /(o = 10. 6/(0 while Eq. (54) takes on the form indicated, 
s 



variation of u vs R_/6 can be generated using Eqs. (65) and 
(70) . An example o¥ this procedure is shown in Figure 6. 

Once the value of tp which best fits the data has been selected, 
the corresponding value of ir can be obtained from Eq. (67) by 
recognizing that R- is determined by the second of the rela- 
tions (69) . Finally, the value of the VP skin friction 
coefficient associated with this profile is obtained by using 
the correspondence (A-9) . 

First Order Closure with Pressure Gradient (FCP) : Since 

V is also assumed to be zero here, the procedure is identical 
to that for ZC in view of the fact that Eqs. (65), (67) and_ 

(70) retain identical forms despite non-zero values of dp/dx. 

First Order CIosutp with Transpiration (FCT) : - As has 
been noted above, the variation of u with R_/cr depends here 
not only on a choice of <p but on a value of^V as well. This 
is true both for Eq. (66) and Eq. (39) . The additional 
relation which is required is Eq. (19) . That is, with W 
specified, is known. Accordingly, with a value of cp 

selected, the appropriate value of V w is obtained from 


V 

w 



(71) 


which follows from Eq. (19) and the definition of Z. With <n 
and V thus selected, the proper value of cr is obtained from 
Eq. (^9) using (50), (51) and (54). Equation (66) is then 

utilized to obtain the desired variation of u vs R-/cr. The 
value of it follows from Eq. (68) while Cj is obtained as 
before. An example of an FCT Clauser Plot is shown in Figure 7. 

Second Order Closure (SC) : The procedure for the SC and 

FCT cases is similar except that the appropriate value of V 
must be determined by iteration using Eq. (36) which, for W 
specified M e and W is of the form H (p,tf w ) = 0>. Otherwise, cal- 
f culation of u vs R -/cr and tt for any value of ~n is identical. 

Figure 8 shows the result of applying SC rules to particular 

■> experimental profiles. A comparison of the results obtained 

using the various methods is presented in the next section. 

A summary of the Clauser Plot procedures outlined above is 
presented in Table II. 
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IV. RESULTS AND DISCUSSION 


A. Profile Analysis 

As indicated in previous sections, one of the primary 
reasons for modifying the compressibility transformation was 
for the purpose of eliminating the observed distortion in the 
wake region of the transformed velocity profiles. This phenomenon, 
which was first pointed out in Reference 3 and also discussed in 
References 2 and 7, manifests itself in values of the wake para- 
meter ir which are substantially smaller than those which are 
expected to occur for constant pressure CP flows. 

To determine whether the modified transformation provides 
any improvement in this respect, a series of profiles were 
analyzed according to the procedures set down in the previous 
sections. That is, Clauser plots for selected profiles were 
prepared and the values of c^ and n determined. A typical set 
of these Clauser plots is shown in Figures 6 and 8 corresponding 
to two choices of closure. The final estimated values of c^, 

7 r, etc., are compared in Table III. As can be noted therein, the 
values of n inferred by using the SC rule are only slightly 
higher than the ones deduced by means of the ZC rules. Figure 9 
shows how both of these values compare with the anticipated CP 
values. It appears that, although some improvement has been 
obtained, the distortion for the most part has not been eliminated. 
This situation prevailed for all of the profiles analyzed which 
included most of those listed in Table II of Reference 2. 


B. Zero Pressure Gradient Results 

To assess the usefulness of the form of transformation 
developed here in predicting boundary layer behavior for the VP 
zero pressure gradient case, a series of calculations were 
carried out corresponding to the range of conditions 0 £ M ^6, 
0.25 ^ W s 1.0. Similar calculations were also carried out 
using both the FCT and FCP forms of the transformation. To 
expedite comparison of these results with experiment, the pre- 
dictions due to Spalding-Chi (Reference 13) were utilized to 
represent the experimental data since this methodology is generally 
recognized as a reasonably accurate correlation of most of the 
available data. 
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The comparison of all of these results is shown in Fig- 
ure 10. From this comparison it would appear that for the 
range of parameters examined: 

. The FCT form of transformation provides the 
best agreement with experiment. 

. The FCP fbrm is superior to the SC form for 
the adiabatic wall case W = 1. 

. For moderate Mach number level (M ^ 2) , the 
SC form of transformation provides improved 
agreement with experiment, relative to the 
FCP form, for flow configurations involving 
heat transfer. 

. For Mach numbers above approximately M =4, 
the FCP and SC forms give essentially the same 
results for the heat transfer case. 


C. Pressure Gradient Results 

The present methodology has been used to generate boundary 
layer predictions for two cases involving adverse pressure 
gradient. These include the adiabatic wall configuration of 
McLafferty and Barber (Reference 14) as well as the case examined 
by Kepler and O'Brien (Reference 15) which involved heat trans- 
fer as well. These results are compared with the experimental 
data in Figures 11 and 12 as well as with the predictions due 
to the FCP form of the transformation. As may be noted, the two 
predictions are essentially equivalent and neither agrees too 
well with the data in the region of maximum pressure gradient. 
Note, however, that this latter disagreement is most probably 
due to the normal pressure gradients which exist in the experi- 
mental setup and which are not accounted for in either of the 
theoretical calculations. 
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V. CONCLUDING REMARKS 


On the basis of the results cited it would appear that the 
modified compressibility transformation developed here provides, 
at best, only marginal improvement over the methodology pre- 
viously developed in Reference 1. This conclusion is valid 
only if the latter is utilized without the use of the so-called 
"Wake Parameter Correlation." In this case the earlier formula- 
tion is superior to the present method. 

In principal, of course, the aforementioned correlation 
could also be employed in the present analysis. The resulting 
formulation, however, can be anticipated to be equivalent to 
that of Reference 1 with the correlation. Thus, it is concluded 
that such additional development would not be appropriate and 
it is recommended that, for the present, the method of Refer- 
ence 1 be employed to treat VP flows with pressure gradient. 

Although the present effort has not been entirely success- 
ful in achieving an improved form of transformation which is 
self-consistent, the basic objective remains valid and further 
effort in this area is recommended. In this connection, the 
approach outlined in Reference 16 appears promising and its 
further exploitation to achieve this objective is indicated. 
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APPENDIX A 


CORRESPONDENCE BETWEEN CP AND VP VARIABLES 


dx . . . 

= « (x) 



(A-l) 

i if - ^ 



(A-2) 

-a(x, 

*-*w 



(A-3) 

U = ~ U 

r? 



(A— 4) 

pv-p% = | t ( P V "Pw V w } + 

ri 5x a 

da-, 

dx J 

(A-5) 

9 = | 

V 



(A-6) 

1 & 1 
6* = ^ [6'* + J (~ -D d y 

] 


<A-7) 

l & l 

H = H + f J (~ -Ddy 

a O P 



(A-8 ) 

C r = P M C r 

f 'w'w f 



(A- 9) 

Ry 

R- = ct J p dRy 

y o 



(A-’IO) 

Rc 

n 0 ^ 

R- = CT J p dRy 



(A- ID 

R e = 5 R e 



(A— 12 ) 
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APPENDIX B 


ELEMENTS OF THE SYSTEM (33)* 


A1 

A2 

A3 

A9 


= lA, 


(i) 


Ai l} + Ai 

+ (2-u*, + 


(i-u*) n* 


2V V r H V <t? -h 

- 1/k i + p- C < 1+ ^ + -f-) 

^ w e e 


A10 
All 

A13 = 


F 

F 

9H 


( 1 ) 

1 

<T?*) 

1 


hlnco 


A14 = 


9-tnH 

9£n<3 


A17 = 2A 


A18 = F 


1 

( 1 ) 


A19 
A2 5 


= F ^*> 


v 2 

- A9 - <p(l + -g— ) 
e 


A26 = 2A 2 (1) + 3A-J 15 - F 2 (1) 


A27 = 2A 2 ( ^ } + (4-u*) A^* } + (1-u*) n* - F 2 ( ^ } 


* The parameters k n , A^' , A^*'* , A^* ) 


( 1 ) . ( 1 ) 


P (D F (r?*) 

2 ' 2 


1' * a ] ' 

(i) _ 0}*) 


u* 


(!) _(rj*) 


are defined in Part II of this report. 
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A28 


1 


A29 

A30 

A31 
A3 3 

A34 

A35 

A3 7 

A38 

A 46 
A47 = 
A48 = 
A55 = 
A56 = 
A61 = 

A62 = 
A63 = 


5H 1 

5-tnU 


a^nH„ 

£ 

a^nU 


= - 1 


: (A1-A25) /V 


w 




F 2 (t? *)/V 


w 


aH, 


av 


w 


a£nH„ 

z 

av 

w 


= - 1 
1 

H + 2 


- 1 

- (H+l) 

aH 

ax" 

a-tnH 2 

ax 

d-tnU /v 
e 

d X 
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A64 


C2 = 


C3 = 


C4 = 


dtnU d-tnu 

(H+l) - (H-H) — - (H+2) 

d x 


d X 


d X 




p w^w CT CT o , J2_ v 

+ ?s ( \ ’ 


V 

1 / i_ + _W ) 

' 5 ( “> S U e 


V u* 


1 f-* _ 1 _ \ 

R V (0 s U J 


G, U ~ 

V 1 e T] 

W w , O a ~ ~ \ 

JT - + { P (p >» a a °’ 


, d£ntJ 
-1 e 
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APPENDIX C 


THE PARTIAL DERIVATIVES APPEARING IN EQS. (36a) and (39a) 
From Eq. (36) 


SH I 2Z 

r ~ ~ % 

b u b u 
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f z 

2 bG 1 

. S 
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» wj 1 
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4* 
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ac 
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« 

<U 

1 
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. e " 
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e 
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where, from Eq. (53), 

as. 


au 

SG 


f ’ \ [' 
= \ [■ 


aT 


2 m 


Y r-f f , 


au s (l-W)+2m e U s j 


aT 


2u 


p u 
B' S 




am^ (l-W)+2m e u 


aT 
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aw 

&G, 


aT 
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where, from Eq. (50) 


ST (l-W)+2m u 
s _ e s 

S vT ~ 1+m 


ST u -T 
s s s 

— TO = TO 

Sm 1+m 

e e 


SW 1+m 


397.2/T 
e 

T 3 - (198 . 6/T ) 2 
s e 


s s 

and where - and =~T — are evaluated in Appendix D. 

hJf.rup S£nv + 


From Eqs. (39) and (54): 
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5-tnHa \ ( 5U s 
Su / lUnv+, 


1-W dXi + ze dX^ 

\ <n d^nv + <a 3 d-f,nv + 
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9X 


5 -t-nHg 
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Where, taking into account also Eq . (51) 
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9m 


9m W+ (1-W) (Xi /<p) +m (X 8 /<p 2 ) 1+m 


9T 
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9W 
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9W W+(l-W) +m (X 2 /<p 3 ) 


9-tnHa _ 198.6 


(T s -1) 
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T -331/T 
s e 


T +198. 6/T 
V S € 


and where “ — 7 and — — =—7 are evaluated m Appendix D. 
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APPENDIX D 


AUXILIARY CP PARAMETERS 


Using Eq. (10) of Part II of this report 
Y + , n / u + 


. 1 f J s -+ - 1 

Xi s -=T J u dy+ = ■=— 

1 y+ v + 

so w 


y S + 


- 1 


(D-l) 


x * s $ s (u+) 8 d?+ = ( 

so w \ 


where 


u + c 


2v 


+ X i 


(D-2) 


V H v> = rr [e x p(v y s +)_1] 


w 


(D-3 ) 


From (D-l) and (D-2) it follows that 
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where, from Eq. (15) of Part II of this report 


d£ny s + [y s + exp (v+ y+ /2)-k a ] 

dlnFT " 1 “ k iy s + exp(v^y s + /2) 


(D-8) 
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APPENDIX E 


LIMITING FORM OF THE ELEMENTS OF THE SYSTEM (33) 


In view of the results shown in Figure 1, it can be antici- 
pated that, in general, V can become vanishingly small during the 
course of any particular calculation. Since several of the 
coefficients appearing in Eq. (33) take on indeterminate forms 
when V =0, special treatment of these is required. Specifically, 
for V W sufficiently small, special limiting forms of these 
coeffYcients are employed by the computer program. 

These special forms are listed below:* 


lim 
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... 3 dv w + 3 s dv w + 
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W+ (1-W) (Xi/^)+m e (X 2 /(0 2 ) 


Limiting values of A34*and A3 5 are not listed here since the 
parameter F 2 and F are themselves proportional to V 
which allows eliminati* — -•*- --- --- 


:ion of the indeterminancy in an obvious 


manner. 
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^ = vii 

dv + 2 \ dv + 3 J 

w w 

V f 2 ^l + 5_ ?t = 

dv + 3 1 d? + 16 Y s 

w w 


dy + y + 3 r . s 

s = Is f _kj_ 1 

d^ w + 4 l l-k 1 y s + J 

Note that these relations and others appearing in the various 

coefficients must be supplemented by limiting values of y g + and 

u . From Eq. (D-3) it is easy to show that 
s ^ 

„ . y + 

lim — lim — _ s 

V -*0 U s V +-*0 U s " (fl 

w w ^ 

while from Eq. (15)^ of Part II of this report it can be shown 

that in the limit V -> 0 y + satisfies the transcendental equation 

w s 

V - ks + k V 


which, for the choice of constants used here (i.e., k a - 4.9, 
1A X = 2.43) leads to 


lim 
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TABLE I. BEHAVIOR OF CLOSURE RULES 
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FIGURE Is EFFECT OF MACH NUMBER, WALL TEMPERATURE AND 
SUBLAYER EDGE VELOCITY RATIO ON THE CP 
BLOWING PARAMETER ACCORDING TO SECOND 
ORDER CLOSURE RULES 





FIGURE 2: EFFECT OF MACH NUMBER, WALL TEMPERATURE 

RATIO AND VISCOSITY EXPONENT ON THE CP 
BLOWING PARAMETER ACCORDING TO SECOND 
ORDER CLOSURE RULES 








FIGURE 4: CLAUSER PLOT FOR CP IMPERMEABLE PROFILE DATA OF REFERENCE 17 





.ocity ratio 








velocity ratio 



FIGURE 7: CLAUSER PLOT ACCORDING TO FIRST ORDER CLOSURE WITH 

TRANSPIRATION FOR PROFILE DATA OF REFERENCE 19 - 
Run 29 (c = .00218, R fl = 6470) 
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FIGURE 10: VARIATION OF SKIN FRICTION COEFFICIENT WITH MOMENTUM THICKNESS REYNOLDS NUMBER, 

MACH NUMBER AND WALL TEMPERATURE RATIO ACCORDING TO THE SEVERAL CLOSURE RULES 
AND COMPARISON WITH EXPERIMENT: SOLID LINE - SPALDING CHI (REFERENCE 13)? 

CIRCLES - FCT PREDICTION: SQUARES - FCP PREDICTION; TRIANGLES - SC PREDICTION 




Factor, H Momentum Thickness, Q -ft . Mach Number 









Form Factor, H Momentum Thickness, 9 -ft. Mach Number 










